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We solve exactly a Terrace-Ledge-Kink (TLK) model describing a crystal surface at a micro- 
scopic level. We show that there is a faceting transition driven either by temperature or by the 
' chemical potential that controls the slope of the surface. In the rough phase we investigate thermal 

fluctuations of the surface using Conformal Field Theory. 
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In this paper, we reexamine a model of crystal surfaces in three dimensions which is inspired by the Terrace- 
Ledge-Kink (TLK) Q ideas of Kossel and Stransky || . By careful specification, an isomorphism can be established 
between our system and an extension of the six- vertex models [Q , which allows an exact discussion at a thermodynamic 
level of the different types of phase transitions which take place. We complement this with recent relevant deductions 
from conformal field theory, which give microscopic insight. We shall return to this six-vertex isomorphism after the 
next few paragraphs, in which we describe the model and outline the results. 

Consider a vicinal section of a crystal surface, that is one which on average is tilted by a small angle from a closed 
packed plane. Let this underlying plane have normal (001) and let the normal to the mean vicinal surface lie in 
the plane containing (001) and (101). Following Kossel and Stransky ||, we give the vicinal surface a microscopic 
O , structure by regarding it as the upper surface of a set of unit cubes, each of which representing a molecule or atom, 
which are stacked vertically with no voids, thereby covering the basal plane (x, y, 0) with columns of cubes. We shall 
| make a somewhat eccentric choice of lattice structure, which is described in detail below, and which will turn out to 
. be of crucial importance in the statistical-mechanical modelling of the short-ranged interactions between ledges. 
' We now impose the vicinal condition by making the height on the left extremum of the surface equal to zero, and 
the height on the right one equal to some constant positive integer. These extrema intersect the lines (0, y) and (M, y) 
of A as zig-zags or zippers, as shown in Fig. ^. Connected components of the upper surface having the same height 
are termed terraces and are separated from other terraces by ledges, which can have bends through an angle of 7r/2: 
these are termed kinks. 

In this paper, we shall make two further restrictions: firstly, there are no adatoms or pits on the terraces; and 
j*^ i secondly, ledges cannot separate adjacent terraces with a height difference greater than unity. This means that 
different ledges can never coincide and that no ledge can form a closed loop: each ledge begins at the bottom of the 
£h , lattice and exits at the top. 

We emphasize that the above construction is not the only way of generating exactly solvable models. Other 
routes are the BCSOS ||U and the surface deconstruction J7| model. There is also the Hamiltonian limit type of 
approximation proposed by Jayaprakash et al. § and by Burkhardt and Schlottmann || , together with other work 
by Villain et al. |l(]] and by Bartelt et al. pi] ]. There is a general review of morphology of vicinal surfaces in p2[ . For 
small tilt angles 9, a polynomial approximation is given to the surface free energy f(0) where 
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f(o) = /(o) + a0 + be 3 (i) 



H ' 

in which a is related to the kink energy and b comes from ledge-ledge interactions which can be entropic in character, 
elastic (although the distinction between these two is not entirely clear) or coming from electrostatic dipolar inter- 
actions in metals. We scrutinize ([!]) from our exact solution and also analyse the ledge-ledge interaction from first 
principles in appendix E. 

Our model shows the formation of a crystal facet at half-filling, which we can regard as at a tilt angle of 7r/4. In 
the low-temperature phase, the surface height fluctuations about the 7r/4 tilted plane are of order one, but above the 
faceting temperature there is power-law decay of ledge-ledge correlation functions, making Kosterlitz-Thouless Jlq] 
behaviour extremely likely. As we shall see, for a certain nonzero repulsion between ledges we get a free-fermionic 
structure of the transfer matrix, and in this case the height fluctuations are shown to be of Kosterlitz-Thouless type. 
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Finally, we note that our model is closely related to a lattice regularization of the anisotropic principal chiral field 




II. THE MODEL 

We consider a square lattice tilted by 45 degrees (see Fig. |l| (a)), and define a configuration by declaring each link 
between two lattice sites to be occupied (state up, denoted 1) or empty (state down, denoted 2). We then impose 
the six-vertex constraint or ice rule: at each vertex (centered at lattice sites) the number of occupied links must be 
conserved between top and bottom. An example of an allowed configuration is shown in Fig. |l](a), where the occupied 
links are drawn as thick lines. 




FIG. 1. (a) six- vertex model configuration; (b) definition of ledge segments for a given six- vertex configuration. 

Every such allowed configuration can now be mapped onto a configuration of terraces, ledges and kinks defining 
a surface. We first introduce ledges by drawing short vertical lines through the middle of each occupied link (see 
Fig. |l|(b)). These segments are then connected by horizontal lines representing kinks, as shown in Fig. |^(a). Treating 
ledges and kinks as steps on a surface, the above procedure then defines terraces, which are indicated by different 
shades of grey in Fig. ||(b). Figure || explicits the three-dimensional structure. 
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FIG. 2. Square lattice configurations in the TLK model: (a) ledges and kinks; (b) terraces. 

The above procedure leads to a number of restrictions on the allowed configurations of terraces, ledges and kinks. 
For example ledges cannot form closed loops, and there cannot be multiple kinks. 




FIG. 3. A configuration in the TLK model with the hexagonal lattice. 

This last restriction can be effectively removed by a slight modification of the ledge representation. Stretching the 
lattice by v3 in the horizontal direction we can introduce a lattice of regular hexagons, and use the sides of these 
hexagons to define the ledges and kinks as shown in Fig. |3|. In this case the restrictions on multiple kinks are built 
into the lattice structure. 




FIG. 4. Section of a crystal surface: (a) square lattice; (b) hexagonal lattice. 

We therefore end up with two slightly different geometrical representations for the configurations: in Fig. |J (a) the 
underlying lattice is square, although there are restrictions in the allowed configurations which involve the tiling of 
Fig. [|; in Fig- 1| (b) the terrace-ledge-kink structure is placed on a hexagonal lattice. 

We now associate a weight with each configuration. 
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FIG. 5. Six-vertex vertices and the two possible interpretations in terms of ledges. 



In terms of the underlying six-vertex model we associate a Boltzmann weight according to the rule shown in Fig. |5| 
with each vertex. The weight w is a kink fugacity, and lower values of w correspond to stiffer ledges. The weights also 
depend on the interaction energy b between neighbouring ledges and the ledge stiffness (note that in our model only 
neighbouring ledges interact and that by construction ledges can never cross each other) . The parameter a is a ledge 
fugacity or chemical potential, which controls the number of ledges, or equivalently the slope of the crystal surface. 
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FIG. 6. The three different assignments of vertex weights discussed in the text. 

We now discuss the relationships between the three weight assignments shown on Fig. |[ Some aspects of the Bethe 
Ansatz for (ii) have been discussed in [ffl. The Bethe Ansatz generates eigenvectors for the (1,1) transfer matrix with 
fixed fraction of up arrows, corresponding to a given vicinal section angle. The solution gives the limiting free energy 
per unit area, in principle for all angles, but in closed form by Fourier methods only for the half-filled case. It was 
argued that there are two quite distinct types of phase transition, which occur at the solutions of e b — (1 ± w) 2 . 

For b > (attraction between ledges), there is a reconstructive phase transition at a temperature T(b) to a state 
with zero limiting entropy for T < T(b). As the temperature approaches T(b) from above, the specific heat diverges 
as (T - T(b)y 1/2 (cf. end of appendix E; this is similar to Q). 

In this paper, we focus on the transition at e b = (1 — w) 2 , for b < 0, implying repulsion between the ledges. This 
problem was discussed in the fixed vertical polarization case with the (1,1) transfer matrix in There it was found 
that at half-filling the system displays an infinite-order phase transition similar to the Kosterlitz-Thouless one Jl5| , 
or the van Beijeren one in the BCSOS case ||. 

We consider the weight assignment (i) of Fig. ^ and we use a Legendre transform with respect to the variable a, with 
no fixing of the polarization, in order to go to the physical situation of case (ii) with fixed polarization. Notice that 
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neither (i) nor (ii) is equivalent to the usual (h,v) field scenario in the six- vertex model (iii): new results are needed. 
We note that the model (i) is not the same as (ii) because the chemical potential a is not always an invertible function 
of the ledge number. In other words, the model [EJ contains phases that are not accessible by tuning a in the model 
considered here. The situation is analogous to the six-vertex model in a vertical field or equivalently the Heisenberg 
XXZ chain in a magnetic field H . If we consider the ferromagnet in the sector with fixed magnetization the ground 



state is rather complicated 1 16 . This phase is not accessible by tuning the magnetic field and considering the absolute 
ground state (i.e. the minimal energy state in the grand canonical ensemble where the number of up/down spins is 
variable) as the ground state will always be the saturated ferromagnetic state with all spins up or down. 

We use the Quantum Inverse Scattering Method, which is reviewed in jlT]], to derive new results, as well as the old 
ones in a formulation which may be more familiar for most readers. We calculate: 

E(a,b) = A lim o ^logZ A (a;&, W ) (2) 

where 

|A| 

Z A (a;b,w) = TraccY[e bni e aini - n2) w n3+ni (3) 
j'=i 

where the trace is taken over all ice-rule-compliant configurations on |A| = MN points for a M x N lattice; rij is the 
number of vertices of type j — 1, 2, . . . , 6 in the lattice configuration. The free energy for fixed polarization is then 
given by 

F(m; b, w) = sup(_E(a, b) — ma). (4) 

a 

Correlation functions may be calculated in either ensemble by appropriately relating m and the optimal a — a(m) in 
(||), which is unique provided that dE/da exists at a — a(m). 

The paper is structured as follows. In section III, we give the Bethe Ansatz equations and expressions for the 
eigenvalues of the transfer matrix. A detailed derivation of these results, as well as the classification of the model 
in the framework of the Quantum Inverse Scattering Method, appears in appendix A. In section IV we discuss the 
phase diagram and calculate explicitly the free energy. In section V we give a precise definition of height and height- 
difference correlation functions. Sections VI and VII are concerned with the determination of the large-distance 
behaviour of these correlation functions. Finally we summarize our results and discuss their relevance in relation to 
other theoretical and experimental work. 

There are several further appendices discussing some technical aspects of our work: the free-fermion case is solved 
explicitly in appendix B, and appendix C deals with the analysis of the Bethe Ansatz equation in the case of negative 
kink energy, where novel peculiarities concerning the distribution of roots in the complex plane are encountered. Ap- 
pendix D gives the second-quantized form of the transfer operator, and in appendix E we discuss effective interactions 
between ledges and some low-angle expansions. 



III. BETHE ANSATZ EQUATIONS 

Denoting the weight matrix by R we have in case (i) : 

R\\ = e a+b , Rll = e~ a , R\l = w, RH = w, R\l = 1 , E*\ = 1 . (5) 

The vertex model defined by (|5|) is exactly solvable by Bethe Ansatz. In a particular parametrization of the Bethe 
Ansatz equations (for the case a = 0) was derived using coordinate space techniques. In appendix A we show that 
the model (|) can actually be embedded into a family of commuting transfer matrices and derive a parametrization of 
the Bethe Ansatz equations in terms of entire functions. In the language of the Quantum Inverse Scattering Method 
p7[ the model (^) is simply a particular case of an inhomogeneous asymmetric 6- vertex model |l^|l|,[n| . As is shown 
in appendix A, it is convenient to introduce parameters 7, lo and v such that 

b _ sin 2 (7-tj/2) _ sinu/2 _ 2v sin( 7 - uj/2) 

sin 7 sin 7 sin 7 

Eigenstates of the transfer matrix (with 2N — M ledges) are parametrized in terms of rapidity variables Afc, which 
are subject to the following set of coupled algebraic equations, called Bethe Ansatz equations 
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( sinh (A fc + i\) sinh (X k - t^) A N sinh (A fc - Aj + ry) , „ , , 

\swh{\ k -q) sinh(X k -i^)J 1 = 1 sinh (A* - A 3 - - i 7 ) ' ' lj 

The eigenvalues of the operator of translation by two sites are given by 

ex P apaxA)) = n ^ fc :1l sin h h ^"SI . («) 

fcVi Smh ( A /c + 1 2 ) Smh ~ 4 2 ) 

and the eigenvalues A,^ of the diagonal-to-diagonal transfer matrix Tdd are 

Add({Aj}He I sin 7 ) lJL s inh(A,-*i)sinh(A,+*^) ' (9) 



IV. PHASE DIAGRAM 

The parametrization (|^) is such that 7 and w cannot always chosen to be real. Furthermore, due to the fact that a 
and b must be real and w must be real and between and 1 there are restrictions on the regions of "physical" 7 and 
lu. As a and b are real so is v, and by exchanging up and down spins in the Bethe Ansatz solution (see appendix A) 
we can restrict ourselves to the case v < 0. 

The permitted values of b and w are covered by the following three regimes of 7 and u>: 

• Region (la): 7 and uj are real with < 7 < ir/2 and —27 < u> < 0; 

• Region (lb): 7 and w are real with ir/2 < 7 < tt and — 2ir + 2j < u> < 0; 

• Region (2): 7 and &> are purely imaginary. It is convenient to choose a parametrization obtained from (^J) and 
(0) by taking 7 — > 17, a; — > ia; and A — > i(A + 7r/2), which leads to the following set of equations 



h 



,, sinh 2 (7-a>/2) sinhtj/2 _ a _ 2l , sinh( 7 - a;/2) 

e = —2 , w = — , e =e — , (10) 

smh 7 smh 7 smh 7 



where 7 > and —27 < uj < 0, 



cos (A fe +2^) cos (A fe _ t~t sin (Afc - A, + fry) fc = 1 M 



COS 



(A fc - if) cos (A fc - »£2±I) y sin ( Afe - Aj - 17) ' 



A dd ( { A,}) = e-(— ) ( ' JJ C0S ^ + ^ C0S ^-^ 
\ smh 7 J ax cos (Aj - z^j cos (Aj + i 11 ^) 



• Region (3): w and 7 — 7r are purely imaginary. A convenient parametrization is then 

b _ sinh 2 (7-cj/2) _ sinhw/2 _ a _ _ 2v sinh(j - w/2) 

smh 7 smh 7 smh 7 



where 27 > uj > 0. We substitute A — * i(X + 7r/2). The Bethe Ansatz equations become 

/sin(Afc+i|) sin (Afc -i 1 ^)^ _ -^t sin (Afc - Aj + £7) _ 
^sin (A fc -q) sin (A fc - i^) / = ^ sin (A fc - A, - z 7 ) ' ~ ' 

The eigenvalues of are 

A fA «- „2a(M-A0 f sinh( 7 - |) \ 2M sin(A 3 - sin (A^ ■ - i^) 

Ad ^i)- e { sinh7 J ll sin ( A j-i2) sin ( A . +i a-) 



(12) 



. M . (14) 



(15) 



By analysing the Bethe Ansatz equations in the various regions we can determine the largest eigenvalues of the 
transfer matrix and thus establish the phase diagram of the model. 
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A. Phase diagram for v = 



As in the case of the symmetric six-vertex model the phase structure of the TLK model for vanishing vertical field 
v = is closely related to the parametrization of the Bethe Ansatz equations. We find that there are three different 
phases corresponding to the three different regimes in 7 and uj discussed above. Phase 1 is a critical (rough) phase 
whereas Phase 2 and Phase 3 correspond to massive phases (smooth surfaces). They are the analogs of the ferro (2) 
and antiferroelectric (3) phases in the symmetric six-vertex model. The resulting phase diagram for v = is shown 
in Fig. @. 




FIG. 7. Phase diagram for zero vertical field v — 0. 

The precise parametrization of the phase boundaries for v = is as follows: the boundary between Region (1) and 
Region (2) corresponds to 7 and uj small with cj/7 finite, and 



e b = (1 + w) 2 with w = . 

2 7 

The boundary between Region (1) and Region (3) is for 7 — tt and uj small with uj / (it — 7) finite 



(16) 



e = (1 — w) 2 with w 



2(tt - 7) 



(17) 



B. Phase diagram for v < 



We now discuss the phase diagram as a function of v for various points in the (7, uj) plane. The situation is 
analogous to the one for the XXZ Heisenberg chain in a magnetic field pQ]. As far as the parametrization of the 



7 



Bethe Ansatz equations is concerned, the (7,0;) plane is divided into Regions (l)-(3) as discussed above. At v = 
the phase structure of the model exactly reflects this division of the (7, to) plane. 

• Let us consider a fixed point Pi — (71, Wi) in Region (1). At v = the model is thus critical. As v is decreased 
the phase boundary between Phase 1 and Phase 2 moves in negative b direction and eventually encompasses Pi : 
the model is then massive and we are in the ferroelectric phase. Pi remains in Phase 2 as v is decreased further. 

• Next we consider a point P2 = (72,^2) in Region (2). For any value of v < the model is in Phase 2. 

• Finally we consider a point P3 = (73,(^3) in Region (3). At v = the model is in the massive Phase 3. As v 
decreases the phase boundary between Phase 1 and Phase 3 is shifted in the negative b direction and eventually 
P3 enters Phase 1. As v is decreased further, P3 finally enters into Phase 2. 

The above picture, together with the precise parametrization of the phase boundaries for v < 0, follows from the 
analysis of the largest eigenvalues of Tdd, which is outlined in the next subsection. The change of the phase boundaries 
under a decrease of v is shown in Fig. ^. The corresponding ledge chemical potential tends to decrease the number of 
ledges in the system and as a result the transition to the state with no ledges (Phase 2) occurs at a smaller value of 
ledge attraction b for fixed w. On the other hand the critical regime extends to stronger ledge repulsion. 




FIG. 8. Phase diagram for v — —0.025 (solid lines) vs v — (dotted lines). 

C. Characterization of the Maximal Eigenvalue 

Let us now turn to the distribution of rapidities A corresponding to the largest eigenvalue of Tdd in the three phases. 
It is convenient to define an "energy" E by 
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E({\ k }) = -]n\A dd ({\ k })\ = 2Na- 



M 

E 

fc=i 



(18) 



where the "bare energy" of a single ledge eo(\ k ) is given explicitly below. The maximal eigenvalue of Tdd then 
corresponds to the minimal energy E. In the following we call the state corresponding to the largest eigenvalue of Tdd 
the "ground state" of the system. 

• Region (1): Let us start with the parametrization for the Bethe Ansatz equations and eigenvalues of the transfer 
matrix of Region (1). The contribution of a single ledge to the energy and to the eigenvalue of the 2-site 
translation operator (H) are 



eo(A) 
JpoW 



-log 


sinh (A - 


f iiji) sinh (A — 




sinh (A - 


- iZ) sinh (A + 




sinh(A + ry/2) 


sinh(A — iui/2 


+ z 7 /2) 



sinh(A - ij/2) sinh(A - iuj/2 - ij/2) 
In the following we restrict ourselves to solutions of the Bethe Ansatz equations for which 



= 1, k = 1 



,M. 



(19) 



(20) 



This means that we do not consider so-called string solutions |2lj,|l^.g2| as they are not important for our present 
purposes: the largest few eigenvalues of Tdd are given by solutions of (Q) that fulfill (20). 

Rapidities fulfilling (20) are distributed on the two lines A = x + iir/2 + iuj/4 and A = x + iuj/4. We find that 
eo(x + in/2 + iui/A) < and £q{x + iuo/A) > for any real x. Therefore for v < the distribution of rapidities in 
the ground state is obtained by filling a Fermi sea of rapidities on the line A = x + in/2 + iuj/4 with |x| < xp. 

The Fermi sea disappears if eo(iu}/4) > and the maximum eigenvalue state then contains no ledges at all. 
Thus for sufficiently large negative v we are in the ferroelectric Phase 2. We exclude this case from the following 
discussion. 

In the thermodynamic limit the ground state is described by a continuous rapidity distribution function p(x), 
which is defined by 



p(x k ) = Urn l/(N(x k+ i ~ x k )) , 



N 



(21) 



where X k — x k + in/2 + iuj/4 are solutions of (]?]). The Bethe Ansatz equations turn into an integral equation for 
p{x) in the thermodynamic limit, which is easily obtained by subtracting the logarithm of (|7]) for k and k + 1 
and then turning sums into integrals 



p(x) 



1 

2^ 



dy K(x-y)p(y) = p {x), 



(22) 



where 



K{x) 
Po(x) 



2 sin 27 



cosh 2x — cos 27 ' 
1 dpo 1 ( sin(7 + w/2) 



sin(7 — oj/2) 



2tt dx it \cosh2x + cos(7 + uj/2) cosh 2a; + cos(7 — w/2) 
The ground state energy density is then given by 



E, 



GS 



dy p{y) £o(y + in/2 + icj/4) 



(23) 



(24) 



where Xf is determined self-consistently by the condition 



9 



5xj 



GS 



= 0. 



(25) 



The density of ledges in the ground state as a function of v is then given by 



D 



dx p(x) 



(26) 



Note that D = 1 corresponds to one ledge per two sites, i.e. half-filling because there are 2N sites in s-direction. 
In the square lattice representation, if we assume that the step height corresponding to a ledge is equal to 
two lattice steps, then the tilt angle of the surface with ledge density D is given by tan0 = D. The angle 
corresponding to half filling is then = tt/4, and we denote by Of = arctan(2) the angle for full filling. 

It is easier to work with the so-called dressed energy e{x) defined by 



1 



X /• 



e(x) - — dy K(x - y) e(y) = e (x + in/2 + iw/4), 

^ J-XF 



where we demand that c(xf) = 0. This fixes xf as a function of v and it can be shown that 



E G s 



dy Po(y) 



(27) 



(28) 



In order to calculate the ground state energy one now simply solves ( [27| ) and ( |28| ) numerically (analytic solutions 
are not possible for generic v < 0). As remarked above we are dealing with a critical phase. In order to facilitate 
the calculation of correlation functions we need to determine the finite-size correction to the ground state energy, 
which is related |^3| to the central charge of the underlying conformal field theory p4| . The leading corrections 
are determined by using the Euler-Maclaurin sum formula when turning sums into integrals Pq,[17t ■ We find 



E GS (N) 



dy £o(y + £tt/ 2 + iw/4) p(y) 



e'(x F ) 



6N 2i:p{x F )' 



(29) 



This can be expressed in terms of the Fermi velocity 

de(x)/dx 



Vf 



dp(x)/dx 



e'(x F ) 
2-Kp(x F ) 



(30) 



where p(x) is the dressed momentum 



p(x) 



dx p(x) p (x) 



(31) 



As a result we have that the central charge of the critical theory is, as expected, c = 1. 

• Region (2): For any v < 0, we have that eo(x + iir/2 + iu>/4) > and eo(^ + iw/4) > 0. As strings do not play 
a role, the ground state is the one with no ledges. 

• Region (3): In this region we need to analyse the Bethe Ansatz equations (^). As we already mentioned we 
encounter three different phases by varying v. For v close to zero the ground state is obtained by filling a Fermi 
sea of rapidities A = x + iuj/4 with x € [—ir/2, tt/2]. The ground state energy is given by ( p8| ) with xf — 7r/2, 
where the dressed energy e{x) fulfills 



i r 



dy 



2 sinh(27) 



cosh(27) — cos(2a; — 2y) 



(32) 



but docs not vanish at the Fermi rapidity 7r/2, i.e. e(±7r/2) < 0. The bare energy and momentum that enter 
(p8|) and (^|) are given by 
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Po{x) 



log 



cosh(7 + ui/2) — cos(2ir) 



4u 



cosh(7 — to/ 2) — cos(2ir) 
1 dpo 1 / sinh(7 + w/2) sinh(7 — w/2) 

27r dx 7r \cosh(7 + w/2) — cos 2x cosh(7 — uj/2) — cos 2x 



(33) 



As increases in magnitude xf stays at tt/2 until, for 
The critical value is characterized by 



larger than some critical value v c , it starts to decrease. 



e(7r/2) 



= . 



(34) 



For values of \v\ slightly larger than v c we are in Phase 1. Energy and ledge density are obtained from the same 
integral equations as before, but now Xf < tt/2. Finally, as \v\ continues to increase it reaches a second critical 
value v' c at which a transition to the ferroelectric Phase 2 occurs and the ground state becomes trivial. 

In Fig. |^ we plot the ledge density as a function v for two points in the (7, lu) plane. These points, denoted by a 
and b, are chosen such that they lie on phase boundaries for v = —0.025 (see Fig. |J). For v — point a is within 
the Phase 1 and the ledge density is 1. As \v\ decreases to v — —0.025 the ledge density decreases until it becomes 
zero: there is a transition to Phase 2. For v < —0.025 we remain in Phase 2 and accordingly D = 0. Point b is in the 
massive Phase 3 for v — and as v decreases from the ledge density remains 1 until at v c = —0.025 we enter the 
Phase 1 and D starts to decrease until it reaches zero at v' c = —0.38. For v < v' c point b belongs to Phase 2. 




FIG. 9. Density as a function of the ledge chemical potential v for two points in the (7, o;)-plane. 



V. HEIGHT CORRELATION FUNCTIONS 



We now define a function h(x, y), which measures the height of the surface at position (cc, y), where the coordinate 
system is such that the transfer direction is parallel to the y axis. Due to the lattice structure, h is constant on 
rectangles of size (1,2) centred on points (x + 1/2, y — 1/2), where x and y are integers with the same parity (see 
Fig. [l^). It is thus sufficient (and indeed most convenient) to measure heights between such points only. 
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FIG. 10. Local TLK configuration (for the square geometry). 



Given the way that we defined the geometry of the surface in terms of ledges and kinks, which are directly related 
to the vertex model, it is clear that height correlation functions are nonlocal quantities. We therefore consider height- 
differences. For any point (x, y) where the two coordinates have the same parity, the occupancy of the link centred 
on point (x,y) is related to the height difference d^(x,y) — h(x + 1/2, y — 1/2) — h(x — 1/2, y + 1/2). This implies 
that d—(x,y) = if the link is empty and d-(x,y) = +1 if it is occupied. If x and y have different parities then we 
define d+(x,y) = h(x + 1/2, y+ 1/2) — h(x — 1/2, y— 1/2). Here also d+(x,y) — +1 if the link at (x,y) is occupied 
and d+{x 1 y) = otherwise. 

The knowledge of the functions d± is sufficient to construct the height function h up to an overall additive constant. 

Since we impose periodic boundary conditions (in height differences), the expectation values of d + and c?_ are 
independent of the position. Moreover, because of invariance under the exchange of left and right, we have (d+) = 
(d_) = D/2, where D is the ledge density D = M/N (the system size of the system is taken to be 2N x 2L). 

The correlation function (d ai (x\, yi)d ao (x , yo)), where cc, = (— l) 1+a! <+i/< for i = 0, 1, depends only on X\ — x and 
Vi ~ Ho because of two-step translational invariance and left-right and up-down symmetries. Any height-difference 
correlation function of the model can therefore be obtained from the following connected correlation function, which 
is defined for all integer values of x\ — xq and y\ — yo: 

C c { Xl - x ,yi - y ) = ({d ai - D/2)(d ao (x , y ) - D/2)) . (35) 

Since we have chosen to work with the two-step transfer matrix T e id, wc need to distinguish between even and odd 
values of y. Our convention is that Tdd connects horizontal lines with odd values of y. 
The correlation function for even values 2y of the vertical separation y\ — yo is given by 

C c (x, 2y) = ^Trace ((1 + a{ - D)T y dd (l + o* x+l - D)r^) . (36) 

For odd separations in the y direction we find that 

C c (x, 2y + 1) = ^Tracc (u^(l + of - D)U 2N rl d {\ + o* - D)t^) . (37) 

The roughness properties of the surface are usually deduced from the behaviour of the height correlation function 

H c (x 1 - x ,y 1 - 2/ ) = ^(ft(a;i,j/i) - h(x ,y ) - {x-i - xq)D/2) 2 ^ . (38) 

This function is related to the connected height-difference correlation function C c by 

C c (x, y) = l [H c (x + 1, y + 1) + H c (x - 1, y - 1) - 2H c (x, y)] 

= l - [H c (x + 1, y - 1) + H c (x - 1, y + 1) - 2H c (x, y)] , (39) 

if x and y have the same parity and 

C c (x, y) = i (H c (x + l,y) + H c (x -1,2/)- H c (x, y + 1)- H c (x, y - 1)) (40) 

otherwise. C c is thus a second lattice derivative of H c . 
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Finally we note that for any state \M, a) in the M-ledge sector the following sum rule holds: 

(41) 



'2N \ 

Y f 0- + <-D))\M,a)=0. 

\x=l / 



This implies that for any y: 

2JV 



Y,C c (x-x o ,y) = 0. (42) 



x=\ 



VI. THE MASSIVE PHASE 



In the massive phase the connected link-link correlation function C c {x,y) decays exponentially in x and y. Exact 
expressions are presently not available although they may in principle be obtained by the method of Jl7| , p6[ . 

It follows that H c (x, y) — E x x — E y y — F (where E x , E y and F are integration constants) also decays exponentially 
as x, y — ► +00. 

Using the sum rule ( [i"2"| ) for xq = N and y = we can show that E x = 0: Inserting (^9|) and (^) into we 
obtain 

H C (N + 1,1)- H c (l -N,l) + H C {~N, 0) - H C (N, 0) = . (43) 

Due to left-right symmetry we have H c (x,y) = H c (—x,y), which implies that H C {N + 1,1) = H C (N — 1, 1). This 
exact relationship is compatible with the large- distance behaviour (|43|) only if E x — 0. 

This argument extends to any asymptotic formula similar to (|43|) provided that all the additional terms have a 
derivative with respect to x which vanishes for large values of x. As we shall see below, this is the case in the massless 
phase as well. 

Returning to the massive phase, if both E x and E y vanish, then the surface is smooth, since H c decays exponentially 
to a fixed value F. The width of the interface is of order y/F. 

However, if E y > 0, then the height correlations in the y direction have fluctuations of order y fy 1 similar to those 
of a one-dimensional string. We anticipate that E y = 0, but we have not been able to prove that this is the case. 

VII. THE MASSLESS PHASE 

We now consider the massless phase (Phase 1). We recall that this phase can be obtained with b and w in Region (1) 
or Region (3), in a certain range of values of v, which we have described above. 

In both cases the maximal eigenvalue corresponds (in the thermodynamic limit) to a filled Fermi sea of Bcthc 
Ansatz rapidities, with a density function given by (p3), where the kernel and the bare density are given by ( p3| ) or 

We use conformal field theory |24| to relate the finite-size spectrum of the transfer matrix to the long-distance 
behaviour of correlation functions. Using the transfer matrix formalism, it is possible to determine the spectrum 
of eigenvalues of r^d for L — > 00, N 3> 1 finite. By conformal covariance, the finite-size spectrum in the toroidal 
geometry of the transfer matrix is related to the power-law decay of correlation functions in the infinite plane p7| . 
In other words, due to conformal symmetry in the critical phase, we can calculate the asymptotic behaviour at large 
distances of correlation functions in the infinite system from the finite-size (in N) spectrum of Tdd- 

Using standard techniques [§5[0 , based on application of the Euler-Maclaurin sum formula to taking the thermo- 
dynamic limit of the Bethe Ansatz equations (Q) , we obtain the following expression for the finite-size spectrum: 

P - P = 2p F d + ^{ M+ ~ M ~ + dAM) (44) 

E - E = ^ ( (^) 2 + (Zd? + M+ + M-\ . (45) 

Here the integers d, AM and M are quantum numbers characterizing the intermediate states in the spectral rep- 
resentation of correlation functions. Given that the largest eigenvalue(s) of the transfer matrix can be thought of 



13 



as (excitations over) a filled Fermi sea of rapidities subject to a Pauli principle pq| , these quantum numbers can be 
interpreted as follows: d is the number of particles backscattered ("2p^-excitations"), AM the overall change in the 
number of particles, and the refer to the creation of particle-hole pairs near the Fermi points ±Pf- 

The other quantities entering (|45| ) are the Fermi velocity vp defined in (|3p| ) and the dressed charge Z = Z(xf), 
which is defined in terms of the integral equation 

Z{x) + ±- r dy K{x - y) Z(y) = 1 , (46) 

27r J-xp 

where the kernel K(x) is defined in ( |23| ) or (|32]). 
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FIG. 11. Dressed charge as a function of ledge chemical potential for two points in the (7, o;)-plane. 




-1 



1 



-1 



In w 



FIG. 12. (a) Dressed charge as a function of ledge interaction for fixed values of ledge stiffness w and chemical potential v 
(b) Dressed charge as a function of \nw for fixed b and v. 



Since the correlation functions in which we are interested contain only operators which conserve the number of 
particles, we can restrict ourselves to AM = 0. We then obtain the following form for the asymptotics of the 
height-difference correlation function 



C c (x,y) 



vpy 2 



B cos(2pi?x) 
(s 5 



v 2 F y 2 ) z2 



(47) 



It is useful to recast ( p7j ) in a slightly different form. We consider a line in the (vf sin9, cos 9) direction (9 = is vertical, 
9 = tt/2 is horizontaljand define a height function hg(r) — h(rvF sin9, r cos 9). In the scaling limit height-differences 
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become derivatives. Thus, we expect the lattice equivalent of h' g = dhg(r)/dr to be cosg+ ^- F sm6 d + + 
which in turn implies the following form for the connected correlation function in radial coordinates 



cos V—vf sin t 



Cefa-n) = (h'einWefa)) = Ag\r 2 - n \~ 2 + Bgcos(2p F v F \r 2 - ri \sm9)r~ 2Z . (48) 

Here Ag and Bg are quadratic forms in cos 6* ± v F s'm9, i.e. affine functions of cos 20 and sin 29. 
Cg{r) is related to the height correlation function 

H c e {r 2 - ri )^ / L(n)-Mm)- j(ri-r 2 )v F smd\ \ (49) 



by 

d 2 H L e {r 

therefore the asymptotic behaviour of H is 



Co(r) = ^ (50) 



rTf , s . , Bg cos(2p F v F rsin9) „ „ ,_„, 

ffg(r)=^lar + , 2 2 . 2fl 2 ^ + ^ + F e (51) 

ip F v F sin 

We believe that Eg = for any value of 0. We have already seen that E n / 2 = 0. 

As — > the amplitude of the ?-~ 2Z term is expected to diverge. We cannot prove this since Bg might be 
proportional to (1 — cos 29), but we have no reason to believe that this is the case. 



VIII. DISCUSSION 

The results of [|J were mainly of a thermodynamic character. In the half-filling case, corresponding to a tilt angle 
of 7r/4, there is a phase transition at e b = (1 — w) 2 for the case of repulsive interactions b < 0. This is of Kosterlitz- 
Thouless, or equivalently, Lieb-F type. For tilt angles < 9 < 9f which differ from 7r/4, there is no therm ody namic 
singularity. For 9 = n/4 and small enough temperature, a Peierls proof along the lines of Brascamp et al. {29| shows 
that the fluctuations in height about the 9 — ir/4 surface are of order one, whereas in the free- Fermi case e b = 1 — w 2 
(see appendix B), there are logarithmic height fluctuations, as in the BCSOS model, which also has a free-Fermion 
point. 

The interesting question is what happens for 9 ^ 7r/4 at a level of correlation functions. We show that for 
9 =/= 7r/4 and no binding of ledges, the excitation spectrum is gapless. This gives a power-law decay of height step- 
step correlation functions by using conformal arguments. In addition to the r~ 2 term which leads to logarithmic 
height fluctuations, we find an oscillatory term which decays algebraicly with a nonuniversal exponent, but whose 
contribution to height fluctuations is subdominant. 

At 9 — 7r/4 and below the transition temperature, the transfer- matrix spectrum has a gap, giving exponential decay 
of correlations of local objects. 

There is another phase transition which occurs for any angle < 9 < 9f at the solution of e b = (1 + w) 2 . This 
has zero limiting entropy per unit basal area below the transition temperature T c , and a specific heat which diverges 
as (T — Tc)^ 1 / 2 from above. In this high-temperature phase, we have gapless excitations as before, and therefore a 
rough surface. This phase transition may be thought of as of pinning-depinning type. The low-temperature phase 
has a complete collapse of the ledges in a close-packed strip with axis in the Euclidean time direction. We have not 
analysed the correlation functions in this case. We must point out that the results in this part of the phase diagram 
are almost certainly critically dependent on the nature of the short-range interactions. 

In appendix E, the interaction between two ledges separated by mean distance I is found to have the same l~ 2 
behaviour as the law which would be obtained from elastic theory or from the interaction between electrostatic 
dipoles located along the ledges, which may be approximated by pairwise interactions between neighbouring edges. 
Our model gives a consistent treatment of that part of the elastic interaction which comes from deformation of the 
surface via entropic interactions. This part tends to zero as the temperature T — > 0. But the dipole term does 
not vanish in the same limit and thus, when present, should dominate the low-temperature behaviour. Thus we 
believe that the model considered by Villain et al. [0 for sections vicinal to Cu(lll) is fundamentally different; our 
model only shows roughening at the unique half-filling angle of tilt, thermodynamically speaking, and the microscopic 
picture strongly supports this view. For less-than-half filling, there is extensive degeneracy because of the form of the 
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short-range interaction in our basic model. One might then ask why the I /I 2 entropic repulsion which comes in at 
the coarse-grained level does not break this degeneracy, for any T > (at zero temperature this interaction vanishes), 
and thereby produce a phase transition. Clearly a better understanding of this is needed. 

The recent work on vicinal surfaces of 4 He crystals |30) suggests that there may be a critical tilt angle of the 
vicinal section, below which neighbouring ledges are rather straight and thus do not overlap, but above which the 
ledges wander sufficiently to come into close contact. In our model, such a critical angle does not appear in the 
thermodynamics. However, in order to address the questions raised by the experiments pCI, it may be necessary to 
study the behaviour of ledges at a microscopic level, and possibly also to take into account defects and boundaries. 
Examples are known pl[ of pair functions near surfaces with weakened bonds where there is a transition behaviour in 
the asymptotics, but, on applying a fluctuation sum formula to get a derivative of the free energy, no thermodynamic 
singularity is seen. Thus caution is advisable in drawing conclusions from the thermodynamics alone. 
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APPENDIX A: A COMMUTING FAMILY OF TRANSFER MATRICES FOR THE TLK MODEL 

In this appendix we construct an embedding of the TLK model into the framework of the Quantum Inverse 
Scattering Method JlTj]. We then show how to recover the coordinate Bethe Ansatz solution of The vertices in 
the diagonal-to-diagonal formulation of the TLK model are given by (||). We define a diagonal-to-diagonal transfer 
matrix as follows 

u 2N = KtMlt ■ ■ ■ *KU* ■ (Ai) 

We note that the square of this transfer matrix commutes with the square of the translation operator tr 

[Ul N ,T 2 R }=0. (A2) 

The transfer matrix Tdd studied in [Q is expressed as 

Tdd = r\ Uj N , (A3) 

and the partition function of the TLK model is given by 

Z = Trace {t m ) L . (A4) 

Correlation functions can also be expressed in terms of Tdd so that the problem of solving the TLK model reduces to 
the one of finding eigenvalues and eigenvectors of Tdd- 

Following |?2[ we can reformulate the problem in terms of an inhomogeneous row-to-row transfer matrix of a 6- vertex 
model in a vertical electric field [jl9| . Starting point is the usual XXZ (symmetric 6- vertex) L-operator jlTj 

sinh(A+i^-a-^) 




U»=\ lSm7 I- (A5) 

Here T is a free parameter and a a are the usual Pauli matrices. Note that there exists a so-called shift point at A 
(at which the L-operator reduces to the permutation operator), i.e. 
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1(1 +<) 



Ln(ii)=[ . • (A6) 



2 



2 (1 



The L-operator of the model with a vertical field is then simply 



£„(A) = exp(|<)MA) exp(|<) . (A7) 



1G 



Consider now an inhomogeneous row-to-row transfer matrix of such an asymmetric 6- vertex model on a lattice with 
2N sites 



T(n) = Trace £ 2 jv(m - i^)^2N-i(^)C 2 N-2{^ - £~)£2JV-3(m) ■ ■ ■ £2(1" - £-)£i(/i) 



Trace exp(| ^ <)T(/x) exp(| ^ <) , 



(A8) 



where T(fi) is the symmetric inhomogeneous 6-vertex transfer matrix constru cted from the L-operators L n (\). Here 
the inhomogeneity f2 is a second free parameter. Choosing fi — i% and using (A6) we find 



Z \ Z Z / Q2iV 02iV 



(^_ 2 (i^-i|: 



«2iV-2P2]V-2 ^ ^ ^ ' "2P2 



(A9) 



where 



C 2 n{i—^—) = 



sin 7 



sin 7 
1 



V 



i(7-f ) 



(A10) 



It follows that the inhomogeneous row-to-row transfer matrix is identical to the diagonal-to-diagonal transfer matrix 
if we make the following identifications 



2,sm(7™f) 



sin 7 



sin 7 



sin(7-f) 
sin 7 



(All) 



Note that 4w = 6 + 2a. Consider now an eigenstate of T(/i) (the transfer matrix of the inhomogeneous, symmetric 
6-vertex model) with eigenvalue v(ix) 



T(j 1 )\K) = v(jjl)\L) . 

Then, because the z-component of total spin is a good quantum number, we have 

T(fx)\A) = e v ^- N ^v(fx)\A) , 



(A12) 



(A13) 



where iV-j-j are the total numbers of up up/down spins in the state |A). This means that we can obtain a complete 
set of eigenstates of T(fi) from a complete set of eigenstates of T(fj,). The eigenvalues of T(/x) are given by (we choose 
the state with all spins down as the reference state) 

^ ^ ^s'mh(fi — — i%) sinh(/x — i%) \ N j-j sinh(/i — Xj + £7) 



. 2 

■ sin 7 



/sinh(/i — i% + iZ) sinh(/i + iZ) N x 



J= i sinh(/i - Aj) 
sinh(/i — Xj — £7) 



n 



V -sin 2 7 / ^ sinh(/i-Aj) 

where M = Nf and where the spectral parameters Xj are solutions of the Bethe Ansatz equations 



(A14) 



/ sinh (Afc + ii) sinh (Aft — 



N 



^sinh(Afc-ii) sinh (X k - £^) 



M 

n 



sinh (Afc — Aj + £7) 
L sinh (Xk - Aj - £7) 



, k = l...M 



(A15) 



We thus find that the eigenvalues of U2N — are given by 



A({Ai» = e a( - M - N > 



sin(7-f) 
sin 7 



M M 

n 



sinh (Aj + i^) 



1 sinh ( A j 



(A16) 
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In order to make contact with the results obtained in Q we need to consider the transfer matrix Tdd = T "jp2N- From 
the standard construction of JlTj it follows that the eigenvalues of the operator of translation by two sites t\ are given 
by 



ex P (i P ({x j} )) = n smh ( Afc ~1? sin M Afc '^? 

K U ' £][ sinh (A fc + q) sinh(A fc -iV) 
This then implies that the eigenvalues Add of r^d are 



(A17) 



AdddX,}) = e MM-iv) 



sin( 7 - |) 
sin 7 



2M A/ 



n 



sinh (Xj + iVj sinh [Xj — i- 1 -^) 



\ sinh (Xj - q) sinh (A 3 - + i 



(A18) 



Equations ( |A15 ) and (A18) are equivalent to (11) and (9) of [jjj if we set a = and substitute the function s(k) of jlj] 
by 



S (fc) 



sinh(A — zw/2 + 27/2) cxp(2A + ij — iu/2) — cxp(iui/2) 



sinh(A - if/2) 



exp(2A) — exp(i7) 



(A19) 



APPENDIX B: FREE-FERMION CASES 



For 7 = 7r/2 the two-particle scattering phase shifts on the right hand side of (0) reduce to —1: the particles are 
free fermions. This corresponds to the following constraint on the weights: e b = 1 — w 2 . The model is physical only 
if e b and w are real and positive . Th is implies that b is negative: the interaction between ledges is repulsive. 



The Bethe Ansatz equations ( A15 ) become 

/ sinh(A fc + ij/2) sinh(A fc -iw/2 + ij/2) \ N ^ ( qj-i (L> 
Vsinh(A fc -i7/2)sinh(A fe -iw/2-i7/2)y 1 ' 



Clearly the rapidities Xk are constrained to lie on the lines 3mA = ilu/A and 5rnA = iuj/4 + «7r/2, which we denote 
by C\ and C2 respectively. The state corresponding to the largest eigenvalue of the transfer matrix Tdd is then obtained 
by filling a Fermi sea of (occupied) rapidities on one of these lines. We distinguish two cases: 

• w coth 2v < — 1 : 

In this case there is a Fermi sea on C2 with Fermi level given by cos!l2Af = — i«coth2u. The line C\ is empty. 
This is a critical phase with power-law correlations. 

• -1 < wcoth2t; < 0: 



In this case the reference state itself (i.e. the state with all spins down) corresponds to the largest eigenvalue: 
both lines Cj are empty. This phase is in general massive i.e. the next-largest eigenvalue of Tdd is separated 
from the largest one by a gap. Accordingly correlations are characterized by an exponential decay. 

Note that for v > the state with all spins up can be used as reference state; the distribution of rapidities is then 
given by the same rule as above but with v changed to —v. (We do not consider this case any further.) 

For any < w < 1 all regimes can be reached by varying a. For v = 0, i.e. a = —6/2, the number of ledges in the 
ground state is M = N. In order to study correlation functions it is convenient to switch to the "coordinate-space" 
notation of H. We define a "particle momentum" k and two functions s±(k) by 

ik _ sinh(A + ij/2) sinh(A - m/2 + ij/2) 
6 ~ sinh(A - 17/2) sinh(A -iuj/2- i-y/2) ' 

s±(k)= w ^ k -^ ± e ifc / 2 y / l- W 2 sin 2 fc/2. (B2) 

We then construct fermionic creation and annihilation operators by means of a Jordan- Wigner transformation on 
the spin variables defining the vertex model 
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In the momentum-space representation there are two "bands" of fermions with momentum k because rapidities on C\ 
and C2 with equal real parts correspond to the same value of k. Fcrmion annihilation and creation operators are of 
the form 



N N-l 



3=1 3=0 
JV AT-1 



Cl(fc) = ^(E e ^'4 + E e* 3 . s±(fc) C + +1 J . (B4) 

The inverse relation between momentum-space and position-space operators is 

c 2j = -=L V e^C k . a c 2j+ i = — L V s a (k)e^C k . a (B5) 



% - V ATf- s _ Q (fc) - SQ (fc) 6 C ^ +1 - V ATf- SQ (fc) - s. a (k) Ck '« ■ [Bb) 

k,a k,ot 

The vacuum state of the fermionic Fock space is defined by Ck,±\Q) = 0. Eigenstates of the transfer matrix are 
obtained by acting with creation operators C\ _j_ on the vacuum 

n 

\ki...k n \a i ...a n ) = l[ct. !a .\Sl) . (B7) 
3=1 

The corresponding eigenvalue of the transfer matrix is 

n 

) =l[e 2a [l + w s a .(k j )][l + w/s a .(k j )} I ai • •• On) • (B8) 

3=1 

In terms of the fc, a variables the maximal eigenvalue of the transfer matrix corresponds to the state |0) defined by 
filling a Fermi sea in the + band between the Fermi points —kp and fcf, i.e. 

10) = II C U k )W ■ (B9) 
lfcl<fc F 

Let us now turn to the calculation of correlation functions, using standard free-fermion methods we can derive integral 
representations for the correlation functions of height differences. However, we shall see that the results are much 
more complicated than for the case of the BCSOS model 



1. Correlators of height-differences in ^-direction 

The height difference in y-direction is given by d y (x, 2y) = h(x, 2y + 1) — h(x, 2y — 1), where x is a half-odd-integer. 
The corresponding operator for x = 3/2 is given by 

2N 
3=3 

where the nonzero elements of A are 

,19 W ,91 W .19 W 2 ,,i W 2 . 

A ll=l 2' A 12 = —i 2 ' A =1 2> A l =-1 2- BU 

1 — ur 1 — w 1 — w 1 — ur 
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The correlation function of the height differences (for even separations) is then given by 

C vy {2x, 2y) = (d y (2x - 1/2, 2y - l)d y (2x - 1/2 + 2s, 2y - 1 + 2y)) 

= Trace (^-^,(3/2)^^(3/2)^^) . (B12) 



We note that this formula is valid only for y > 0. formula The operator (BIO) can be expressed in terms of the 
fcrmion operators as 

6 y (3/2) = i ^ a aijaa (fci,fc) Ct a (fe) C ai (fc) , (B13) 



AT 



where 



., , w s ai (ki) + we lkl lk2 s a2 (k 2 ) + e lkl + s ai (ki)s a2 (k 2 ) ,_ >1 

Glai,a 2 (fcl, «2j = 7 7T7 T^T n — v\ • (B14) 

(w - w- x )(s a2 (fc 2 ) - S- a2 (k 2 )) 
This yields the following representation for the correlation function of height-differences 

C vy (2x,2y) = ^ a+Akh,kp) a a ,+(k p ,k h ) {A a (k p )/A + (k h )) y e lx{k -~ kh \ (B15) 

holc(fc^ , + ),particle(fc p ,a) 

Here the sum extends over all holes with momentum kh in the Fermi sea in the + band (\kh\ < kp) and all particles 
in either the + band (tt > \k p \ > kp, a = +) or the — band < tt, a — — ) and 

Aa(fc) = e 2a [l + to s a (k)][l + w/s a (k)} . (B16) 

In the thermodynamic limit we obtain the following integral representation 



C yv (2x,2y) 











dk h < 

















wcos(k p /2) + yl — w 2 sin 2 (fc p /2) 
;cos(fc h /2) + yjl-w 2 sin 2 (k h /2) / 



■2y 

ix(k p -k h ) 



+ I dk h / dk p a + _(k h ,k p ) 

— 7T 



-fc« 



wcos(k p /2) - yjl - w 2 sin 2 (fe p /2) 
K wcos(kh/2) — \Jl — w 2 sin 2 (fcft/2) ( 



2» 

e ix(k p -k h ) ^ (B17) 



We have not managed to greatly simplify (B17). However, it is straightforward to extract the asymptotic behaviour 



for z — > 00 by expanding the integrand around the maxima at kh, k p — ±kp. 

2. Correlators of height-differences in x-direction 

The height difference in x-direction d x (l, 2y — 1) = h(3/2, 2y — 1) — h(l/2, 2y — 1) is measured by the operator 

w-^-ji £ ( , M ,, .V W ' !>M ' (B18) 

The corresponding correlation function of height differences (for even separations) has the following integral repre- 
sentation 
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C xx (2x,2y) 



dk h 



s+(k h ) 



s + (k p ) 



s + (k p ) - s_(k p ) s + (k h ) ~ s^(k h ) 

ix(k p -k h ) 



w cos(fc p /2) + ^Jl — w 2 sin 2 (fc p /2) 
wcos{k h /2) + yT -w 2 sin 2 (fc, l /2) / 

' kF dkh r dk s+{kh) s - {k ^ 

-k F J-n P s-(k p ) - s+(k p ) s+(k h ) - s-(k h ) 



w cos(fc p /2) — \Jl — w 2 sin 2 (fc p /2) 



2;/ 



Jx(k p -k h ) 



w cos(k h /2) — J 1 — w 2 sin (fc^/2) 



(B19) 



3. Height-correlations along the {/-direction 

Above we have derived integral representations for the (local) correlation functions of height-differences. In the 
remainder of this appendix we consider height-correlations along the transfer direction. The height correlation function 
H(Q, y) = (0(3/2, 2y) - h(3/2, 0)) 2 ) is related to C yy (Q, y) by 

H(0, 2y + 2) + H(0, 2y - 2) - 2H(0, 2y) = 2C yy {Q, 2y), (B20) 



where #(0, 0) = and H(0, 2) = C yy (0, 0). Summing (|B20| ) we find that 

y-l 

C yy (0, 0) + 2 ^ C yy (Q, 2z) = H(0, 2y) - H (0 , 2y - 2). (B21) 

z=l 

The algebraic decay C yy (0, 2y) oc ^ of the height-difference correlation functions implies that the surface is rough 

H(0,2y) (xF + E y y + A 1 hxy + ... , (B22) 

where in the thermodynamic limit 

OO 

E y = C yy (0,0) + 2j2C yy (0,2y) . (B23) 

y=i 

On physical grounds we expect that E y — as we now argue. Let us define a variable D y = X^j/=i ^(^i 2y) in our 
finite 2N x 2L-toroidal geometry. is independent of x because d x (x + 1/2, y) = h(x + l,y) — h(x, y) is subject to 
periodic boundary conditions. E y and D y are related by 

(D 2 ) 

E y =\im (d y (x,y)D y )= lim (B24) 

The average tilt of ledges in the xy-plane is s = (N D y ) / (LM): for instance if Z?j, = the ledges are on average 
parallel to the transfer direction (s = 0). Phenomenologically we expect that ledges have a stiffness E, so that the 
total energy cost of a nonzero D y , i.e. of a nonzero s is roughly T,LMs 2 /2, at least for s< 1. The equipartition 
theorem then gives that the average value of s 2 should be of order (NM)^ 1 . This implies that E y = provided that 
in thermodynamic limit iV 2 increases faster than the number of ledges M . This condition is satisfied whenever the 
ground state is described by a finite density function, with M/N finite. 

For the special case at hand we have explicitly verified that E y = 0: Inserting ( |B15| ) together with 

C yy (0,0) = (0\(w- w' 1 ) {wcfc 1 c 2 4+wc 1 ctc+c 2 + c 1 c+ -c+c 2 } |0). (B25) 



into (B23) we obtain 
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+ D a + , a (fc /t ,fc P K, + (fe P ,M A (fc ) _ A (M , (B26) 

holo(fc h ,+),particle(fc p ,a) +V ' K PJ 

where we have explicitly displayed the dependence on the Fermi momentum kp. By inspection we see that 

E y (0) = . (B27) 

Rather than trying to evaluate ( B26 ) directly, we calculate the change NSE y (kp) = N(E y (kp + 5kp) — E y (kp)) under 
a change of Fermi momentum 5k f = 2irN~ 1 

N8E y (kp) = — — 

10 1 — W S + y p — S-,_F 

+ 1Tt («+,+(fcF, fc)a+,+ (fc, k F ) + a +>+ (-k F , k)a +>+ (k, -k F )) ^ + ^ F \ + ^ + ^\ 

+ ^ E («+,-(fcF, k)a-+(k, k F ) + a+,-(-k F , fc)a_ + (fc, -fc F )) ^±|^j ± ^ . (B28) 

Turning the sums into integrals and integrating numerically we find th at th is is indeed zero. Alternatively, if we 
perform the sums for large finite N we find that the right hand side of ( B2§| ) scales like l/N. In conjunction with 
(B27) this shows that E y is indeed zero. 



APPENDIX C: THE CASE OF NEGATIVE KINK ENERGY 



In this appendix we discuss the case w > 1. 

This corresponds to the following extensions of the regions discussed earlier: Region (la) with — 2-7T < uj < —2j; 
Region (lb) with — 2ir < u> < —2tt + 2^\ Region (2) with to < —2j; Region (3) with uo > 2-f. In all cases the one-particle 
spectrum is quite different from the results obtained for w < 1. In addition to the lines in spectral parameter space 
discussed above, we find that there are curves on which the bare momentum is real and the bare energy has a constant 
real part but acquires a variable imaginary part (whereas both bare momentum and bare energy remain real on the 
straight lines). 

We specialize the remainder of this discussion to Region (3). The curves, shown on Fig. [l^, are parametrized by 
A = x + iy + iuj /A with x and y real and satisfying: 

cosh7COsh(2y) = cosh(w/2) cos(2a;). (CI) 
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FIG. 13. Location of one-strings in the (x, y) plane with A = x + iy + iuj/4, for 7 = 1 and uj = 5. 



The line corresponds to low- momentum states — p c < po < p c , and the curves to high-momentum states \po\ > p c 
(p c is comprised between and tt and depends on the parameters 7 and lo). If we consider the dispersion relation, i.e. 
the (bare) energy as a function of (bare) momentum, there are branch points which correspond to the contact point 
between the line and the curves in term of the spectral parameter. For the high-momentum states the bare energy 
becomes complex, and its real part is constant on the whole, equal to the energy at the contact point. 

We first describe the maximum-eigenvalue state. We believe that this state is made up of particles with spectral 
parameters \ k — x k + iw/A, where all x k are real. However the situation is different from that considered in the main 
text in one important aspect: for uj > 27 we find that the bare momentum is not a monotonous function of x, so that 
the bare density po(x) — ^dp /dx is not positive: it turns out that its average value vanishes. 

The logarithm of the Bethe Ansatz equations is: 

Np (x k )-J2 d ( x k-Xj) =n + 2nl k (C2) 

3 

where 

e"<-> = gj£ + *?) (C3) 
smh(ai — ij) 

and the I k are integers, which depend on the cut structure used for 9. We make the following choice: 9(x) is continuous 
throughout the interval — 7r < x < n, with 9(— tt) = 2tt, 9(0) = 0, 9(n) = —2ir. 

A direct numerical solution of these equations for small systems suggests that for this choice of cut structure, and 
provided the spectral parameters are ordered (x k < x k + l), then for the state corresponding to the largest eigenvalue 
of the transfer matrix we have Iu+x — Jfc = 1. 

If we assume that this structure is correct, and take the thermodynamic limit, we find that the density of particles 
p(xk) = (N(xk+i — Xk))^ 1 satisfies the following integral equation: 

1 f XF 

p(x) + K(x- y)p{y)dy = p (x), (C4) 

where the kernel is given by: 

= f- - : 2 o Sinh27 o ■ (C5) 
ax cosh 27 — cos 2x 

Since the bare energy is negative for any value of x we assume that the Fermi level is xf = 7r/2 for v = 0. We can 
therefore solve the integral equation by using Fourier decompositions. The result for the density is then: 

o(x) - 1 sinh(o;/2 - 7) 

P{ ' ix cosh(u;/2 - 7) - cos(2x) K ' 

We have performed further numerical calculations in order to confirm this result. 

First, for systems up to size 16 (i.e. N = 8) we have calculated directly the largest eigenvalue of the transfer matrix 
using an iteration method. This gives a good approximation of the ground state energy per site, and the result agrees 
with the analytical expression obtained above. 

For much larger systems (several hundred sites) we have solved the Bethe Ansatz equations numerically, using the 
distribution of spectral parameters predicted by our theoretical result for the density as "seed" for a Broyden algorithm. 
The agreement obtained by this technique confirms that the integral equation we solve does indeed correspond to the 
thermodynamic limit of the Bethe Ansatz equations. 

Let us comment briefly on "excited states" , i.e. eigenvectors of the transfer matrix corresponding to subleading 
eigenvalues. One type of excitation corresponds to removing a particle from the Fermi sea without introducing a 
momentum. Its energy vanishes for v — —7/2. 

Other excitations can be constructed by creating particlc-holc pairs. However, any eigenstate of the transfer matrix 
must have a real momentum. This forces us to introduce two particles at spectral parameters X p — x p ± iy p + iui/A 



with (x p ,±y p ) on the curve (CI). The two holes have spectral parameters \ = Xi + ilu/A with Xi real for i = 1,2. 



The momentum and energy of this excited state is found to be 
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f(x 1 )f(x 2 ) 

-SB \sm{xi + iw/4 - i 7 /2) sin(x 2 + ius/4 - ij/2)\ 2 
e = T^ - ? T~ T~ — 71 ■ m\ ■ i T - : — Ta i ■ /oM2 ' \^°) 



where 



J v ; sin(x - ?;w/4 + i 7 /2) v ; 



APPENDIX D: SECOND-QUANTIZED FORM OF TRANSFER OPERATOR 

The transfer operator which we use is a product of noncommuting terms, which operate on alternating sublattices. 
Consider a vertex as shown in Fig. |^ which separates two vertically-consecutive horizontal rows of edges. When an 
edge is occupied by a ledge, let it be occupied by a particle, or an up spin. Then the transfer operator for site j is: 

Tj = (1 + e- T (afaT_ t + <7+ +1 aJ)) (l + (e b - l) nj n j+1 ) (Dl) 

with rij = 2erJ — 1. Since n 2 = rij, basis states with ct|ct| +1 = 1 are annihilated by (rfcj+i + <jJ j+i <J J)^ but on those 
with (T|cr|_|_ 1 = —1 we find that ((Tj~crJ +1 + < J ^J r i< J j) 2 acts as the identity operator. It follows that 

Tj = (sinh 2r) 1/2 exp (r* {(rj(Tj +l + exp (b nj n J+1 ) (D2) 

which may be written in terms of Fermi operators in the usual way. This allows us to make contact with the appendix 
of |sj. Note however that, contrary to the statement implied there, Tj and Tj + ± do not commute, because the 
"hopping" terms do not allow it. But, if we take b and t* so small that only linear terms can be retained, where r* 
is defined by w = tanh r* , then the transfer matrix 

N N 



i 

can be approximated by T sa (2 sinh2r) Ar e ff with 

2N 

H = bJ2 njn j+ i + t* [tfaj+i + <?J ) • (D4) 



2N 2JV 



1 



This is the XXZ model in a field. We emphasize that (D4) is an uncontrolled approximation which is not needed in 
our solution. 

APPENDIX E: LEDGE-LEDGE INTERACTIONS 

The interaction between two ledges may be investigated by considering the incremental free energy per unit cylinder 
length with 2p ledges (the 2p-particle sector, that is) on a cylinder having finite circumference 2N, length L and joined 
along the cylinder axis to form a torus. We define the incremental free energy by 

'"<»•»>- &£.Z*%$$ < E1 > 

where Z2 P (L,2N) is the partition function o n the torus for 2p particles. This is evaluated by giving the momenta 



their lowest values (mod 2ir) consistent with (A15). Asymptotically for large N, in that part of the parameter space 
with no ledge binding, this equation becomes 

for j = 1, . . . 2p with minimal solutions 
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From @, we have 



(2j-2p-l)ir/N. 



V a Po J p =0 i 



(E3) 



(E4) 



We assume that the second term above is made up of 2p equal ledge-ledge pair interactions and that the ledges are 
on average equally spaced at a distance I — N/p. Then this asymptotic pair interaction u(p, I) is given by: 



u{p, I) 



7T 2 fd 2 e 



2pH 2 \ dp 2 

TT 2 W 2p 2 - 1 
121 2 p 2 



) X> - 1) 2 

/ p =0 1 



(E5) 



Notice firstly that u(p, I) is independent of b. It depends on the sector p; this is an example of the effect mentioned 
by Fisher p4l in the wetting context. Much interest attaches to the case with a small density D of ledges, which is 
the same as that given by (E5) in the p — > 00 limit. In this regime the following approximation can be used: 



where 



and 



D = 2x F p(0) + P W (0)^ + 0(4) 



P (o) = (i-k(o)^) { 1] Po (o) + 0(4) 

p< 2 )(0) = pW(0) + (^) V)(0) P0 (0) + 0(4) 



V 7T / 



It is natural to define tjf by 



y F = x F [l- K(0) 



so that ufPo(0) = %fp(0), a dressing of the Fermi level. Thus 



D = 2p (0)y F 



P (2) (0) 



Vf + 0{y%) 



which is readily inverted, giving 



Vf{D) = 



1 D p (2) (Q) 

2p (0) 48p (0) 4 



D 3 + 0(D 4 



(E6) 



(E7) 



(E8) 



(E9) 



(E10) 



(Ell) 



The free energy is obtained from ( J24| ) by approximating the small integral for small xf and then inserting (|ES|), 
leading to 



E = 2c 



3 

a + 2y F po(0)e (0) + ^ (e^ 2) po(0) + e p^(0)) + 0(y 



2a + e a {0)D- 



24 p o (0) 2 



0(L> 4 ) 



Now 



4 2) (Q) 

Po(0) 2 



2n 2 w 



(E12) 



(E13) 



giving 
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E = 2a + e (Q)D+—^- + O{D i ). (E14) 

A similar expansion can be used near the transition between the massless and ferromagnetic phases. This transition 
can be driven either by the chemical potential or by temperature. In both cases the analytic behaviour near the 
transition can be studied using the parameter t = eo(0), which is equal to a — a c if the relevant parameter is the 
chemical potential a, and t oc T — T c if using the temperature. The Fermi level is obtained by expanding the left-hand 
side of e{xp) = 0: to leading order e(0) = eo(0) and e^ 2 '(0) = 6^(0), therefore xf oc y/i. It follows that the free 
energy has a i 3 / 2 singularity, and thus that the specific heat diverges like i" 1 / 2 . 
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